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Abstract 

We study the fermion mass and mixing hierarchy problems within the context of the SU{5) 
4+ld domain-wall brane model of Davies, George and Volkas. In this model, the ordinary fermion 
mass relations of SU (5) grand unified theories are avoided since the masses are proportional to 
overlap integrals of the profiles of the electroweak Higgs and the chiral components of each fermion, 
which are split into different 3+ Id hyperplanes according to their hypercharges. We show that 
the fermion mass hierarchy without electroweak mixing can be generated naturally from these 
splittings, that generation of the CKM matrix looks promising, and that the Cabibbo angle along 
with the mass hierarchy can be generated for the case of Majorana neutrinos from a more modest 
hierarchy of parameters. We also show that under some assumptions made on the parameter space, 
the generation of realistic lepton mixing angles is not possible without fine-tuning, which argues 
for a flavour symmetry to enforce the required relations. 
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I. INTRODUCTION 



In the standard model (SM), three of the most open problems are how the fermion 
mass hierarchy is generated, the origin of small mixing angles in the Cabibbo-Kobayashi- 
Maskawa (CKM) matrix and near tribimaximal mixing in the lepton sector. With neutrino 
masses now known to be nonzero but under 1 eV, the mass hierarchy has a spread of at 
least 14 orders of magnitude, given that the top quark has a mass of roughly 170 GeV. 
Amongst approaches used for solving these problems are grand unified theories (GUTs), 
higher dimensional operators, and flavor symmetries. 

One of the most promising new theoretical frameworks for solving hierarchy problems that 
has emerged over the last decade has been extra-dimensional models, such as the Arkani- 
Hamed-Dimopoulos-Dvali (ADD) model [1], and the two Randall-Sundrum (RS) models 
[21 [3]. The ADD and RSI frameworks solve the hierarchy problem between the Planck scale 
and the electroweak scale, which is of a similar order of magnitude to that of the fermion 
mass spectra. For other papers on extra-dimensional models, see Refs. ^-E]- 

In RS2 models the gauge hierarchy problem is not solved by extra-dimensional physics, 
but the split fermion idea of Arkani-Hamed and Schmaltz [9J can be used to generate fermion 
mass hierarchies from exponentially sensitive overlap integrals of extra-dimensional profile 
functions. Similarly, the RSI setup can address this problem by allowing fermions to prop- 
agate in the bulk and thus acquire non-trivial profiles [TUl [TT]. The idea is that the 3+ld 
fermion zero modes are in general localized around different locations along the extra di- 
mension, with dimensional reduction then producing an effective 3+ld Yukawa coupling 
constant that is the product of the 4+ld Yukawa coupling constant and an overlap integral 
involving profile functions. When the profiles are split, the overlap integrals are suppressed, 
leading to small 3+ld effective Yukawa coupling constants. This fits in well with the fact 
that quark and lepton masses, except for the top quark, are suppressed with respect to the 
electroweak scale. Scalar bosons will also in general be split, a phenomenon we shall use to 
suppress colored-Higgs-induced proton decay (see Refs. [T21 - I2T] for more on the use of the 
splitting of fermions and bosons in extra dimensions to generate fermion mass textures and 
to suppress proton decay and other baryon number violating processes). 

In this paper, we shall utilize the SU (5) 4+ld domain- wall brane model devised by Davies, 
George and Volkas (DGV) [22j to address the fermion mass and mixing angle problems. In 
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this RS2-like model, the spht fermion idea arises naturally, and thus the usual SU{5) quark- 
lepton mass relations are avoided. It will be shown that the mass hierarchy problem can 
be solved using this method, and that the mass hierarchy and the Cabibbo angle can be 
accounted for in the two-generation case with Majorana neutrinos. We also explain why 
tribimaximal mixing cannot be accounted for without fine-tuning, and that the addition 
of a flavor symmetry therefore seems necessary. We are thus led to the view that extra 
dimensions provide an excellent way to qualitatively understand mass hierarchies, but they 
are insufficient to explain all the observed mixing angle patterns. The reason the flavor 
problem has proven to be so difficult may be because more than one ingredient is necessary: 
extra dimensions on their own, and flavor symmetry on its own, are only partially successful. 

The following section reviews the DGV model and develops it further in several important 
ways: neutrino mass generation is examined and the see-saw mechanism implemented, and 



the dynamics of scalar-field localization is shown to be analytically tractable. Section III 



then analyses the parameter space of the model to produce the required mass and mixing 
angle hierarchies, with the aforementioned caveat for tribimaximal lepton mixing. Section 
is our conclusion. 



II. THE MODEL 



The DGV model is a specific extra-dimensional theory featuring the brane as a topological 
defect: a kink-like domain-wall configuration [23j. Domain walls are stable classical solu- 
tions of suitable scalar field theories that exhibit a brane-like character, with energy-density 
peaked around the centre of the wall. Unlike fundamental branes, they have a finite width, 
and are most naturally used to replace the 5-function-like fundamental brane of the original 
RS2 model. Like RS2, a 3+ld graviton zero mode is dynamically localized. Unlike the 
original RS2 setup, all other degrees of freedom (fermions, scalars and gauge bosons) must 
be dynamically localized. In contrast to the fundament al-brane case, it is not possible to 
simply postulate that various fields are confined to a domain-wall: one must have dynamics 
to do it, and that is the main challenge in developing realistic models of this kind. 

The dynamical localization of chiral fermion zero- modes is automatic when 4+ Id fermions 
Yukawa-couple to a background scalar field in the form of a kink [23] . Thus the chiral fermion 
structure of the SM can be naturally accommodated. Similarly, additional scalar bosons 
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such as a Higgs doublet can be dynamically localized to a domain wall through a Higgs 
potential that couples those extra scalars to the background scalar field configuration 
Such localized scalars can even obtain negative squared masses, thus triggering spontaneous 
symmetry breaking on the wall. 

The most difficult issue is the localization of gauge bosons, with the need to maintain exact 
3+ Id gauge invariance to ensure gauge universality. A promising mechanism was proposed 
by Dvali and Shifman [25J, the physics of which is quite different from the localization of 
fermions and scalars. The idea is that a gauge group G spontaneously breaks to a subgroup 
H inside the wall, but is restored in the bulk. The bulk gauge theory is taken to be in 
conffnement phase. The proposition is that the gauge bosons of H are then dynamically 
localized to the wall as exactly massless states enjoying exact 3+ Id gauge invariance. There 
are two heuristic arguments for why this should be the case. Dvali and Shifman themselves 
argued as follows: Take the case where G = SU{2) and H = f/(l), and call the gauge 
boson of f/(l) the "photon". The photon is obviously free to propagate as a massless gauge 
boson in the plane of the wall. But in propagating transverse to the wall, into the bulk, the 
SU (2) confinement regime is encountered, and the propagating states must be colorless and, 
importantly, massive glueballs. The photon must incorporate itself into a massive glueball 
to enter the bulk. But this mass gap makes this transition energetically disfavored, thus 
trapping the photon on the wall. Subsequently, Arkani-Hamed and Schmaltz |26] presented 
another heuristic picture: the photon field lines must be repelled from the bulk, because a 
confinement-phase region is by definition unable to support diverging electric fields. Thus 
the fiux is channeled along the wall, effecting a dimensional reduction. At large distances 
within the wall away from the source, the field lines exhibit 3+ld Coulomb form. If the 
source is instead located in the bulk, then a fiux tube leading to the wall is formed, with the 
field lines then diverging outward as if the source was located within the wall. Thus, the 
long-distance behavior of the field lines within the wall is independent of where the source is 
placed. If the source is smeared out along the extra dimension due to a profile function, then 
the corollary is that the asymptotic field line behavior is independent of the profile. Charge 
universality is thus maintained, no matter how the source is distributed along the extra 
dimension. These conclusions generalize to an arbitrary G and H, provided any glueballs 
associated with H (e.g. QCD glueballs) are less massive than G glueballs. 

So, the way to develop potentially realistic domain-wall brane models is clear. One first 
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postulates a scalar field theory that admits a suitable topological domain wall solution. 
This configuration must break gauge group G to H , and H must contain or be the SM 
gauge group. Chiral fermion zero modes and additional scalars such as an electroweak Higgs 
doublet are then dynamically localized as sketched above. The DGV model uses the minimal 
gauge structure where G = SU{5) and H = SU{3) x SU{2) x f/(l) [22]. 

A. The field content 

The scalar fields in the model are 

rj^l, (1) 
X ~ 24, (2) 
$ ~ 5*. (3) 

The fermion content of the theory consists of the SM fermions with a gauge singlet right- 
handed neutrino for each generation i = 1, 2, 3. The SM fermions are placed into the 
following SU (5) representations, 

% ~ 5*, (4) 
~ 10, (5) 

while the right handed neutrinos are singlets, 

N' ~ 1. (6) 

contains the charge conjugate of the right-chiral down-type quark and the left-chiral 
lepton doublet, and "^{^ contains the left-chiral quark doublet and the charge conjugates of 
the right-chiral up-type quark and electron-type lepton for the generation i. 

In this 4 + Id model, matter is confined to a domain-wall brane formed from a solitonic 
kink configuration for the rj field. To implement the Dvali-Shifman mechanism, SU{5) is 
broken inside the domain wall by the second background field x which transforms under the 
adjoint representation. It attains a non-zero value for the hypercharge generator component 
inside the domain wall, so the gauge group respected on the domain wall is that of the SM. 

Chiral fermion zero modes are trapped on the domain wall through Yukawa interactions 
with the Tj and x background configurations. Similarly, additional scalar fields are trapped 
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by introducing quartic interactions between those scalar fields and the background domain 
wall. 

For the purposes of this paper, we shall ignore gravity, although a similar analysis will 
have to be done with its inclusion in a later paper. It has already been noted that the RS2 
graviton localization mechanism also works for a domain-wall brane. For a discussion of 
how gravity affects the dynamical localization of other fields, see Refs. [221 EZ]- 

B. The background domain wall configuration 

The background domain wall configuration is formed from a self-consistent classical so- 
lution for the coupled fields t] and x- The singlet scalar field rj forms the kink-like domain 
wall, while the adjoint x breaks SU{5) down to SU{3)c x SU{2)l x U{\)y on the domain 
wall by attaining a bump-like configuration. 

The relevant part of the action for describing the dynamics of the background is [22], 



where T contains all the SU{b) gauge-covariant kinetic terms for all the fields. V"^^ is the 
part of Higgs potential containing the quartic potentials for r] and x, with 



= (cr?' - ^^l)Tr{x^) + ar^Trix") + \i[Tr{x^)f + \2Tr{x') + /(r^^ - v'^f. (8) 



We want x to break SU (5) to the SM on the domain wall, while having the bulk respect 
the original gauge symmetry. We do this in the standard way by giving the component Xi 
associated with the hypercharge generator Y a non zero value on the brane, and having all 
the other components vanish. Thus the potential reduces to 



where A = Ai + 

To find the background configuration, we need to solve the Euler-Lagrange equations for 
1] and xi subject to the boundary conditions 




(7) 



(9) 



ri{y = ±oo) = ±v, 



Xiiy = ±oo) = 0, 



(10) 
(11) 
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which are degenerate global minima of Vrj^. For the sake of simplicity, we choose to impose 
the constraints 

2fil{c - ~X) + {2c~X - 4fX - c^y = 0, (12) 
a = 0, (13) 

yielding the analytic solutions, 

r]{y) = vtanh{ky), ^^^^ 
Xiiy) = Asech{ky), 

where k"^ = cv^ — u^, and = - '^^ . We should stress that the above conditions are 
not fine tuning conditions, and they are chosen simply so that the background fields obtain 
analytic forms. To find solutions, these conditions need not be imposed, and for a finite 
range of parameters we can always find numerical solutions which are kink-like for rj and 
lump-like for x i22j • The graphs of these solutions for i] and xi shown in Figures [T]^a) 
and[T](b) respectively. 




3" 



(a) (b) 
FIG. 1: (a) The kink r] and (b) Xi functions of the extradimensional coordinate y 

The kink-like ri{y) has its energy density localized about y = 0, forming the domain wall 
brane. The bump-like Xiiu) breaks SU{5) to SU{3)c x SU{2)i x f/(l)y on the domain wall. 
To preserve the topological stability of the domain wall, a spontaneously broken Z2 
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reflection symmetry must be introduced. Under this discrete symmetry transformation, 



y ^ -y, 

7] ^ -7], 

- X, (15) 

5 AT-i 



wliere tlie 4+ld Gamma matrices F*^ for M = 0, 1, 2, 3, 5 are defined 

T'' = Y, forM = /i = 0,l,2,3, 



(16) 



For tlie entirety of this paper, 4+ld Lorentz indices will be denoted with upper case Roman 
indices, while ordinary 3+ Id Lorentz indices will be denoted with lower case Greek letters 
as usual. Also, = y. The next step is to localize the fermions to this background. 



C. Localizing the charged fermions and the left-chiral neutrino 

To localize the charged fermions and the left-chiral neutrino to the domain wall, we need 
to couple them to the background. The relevant Yukawa Lagrangian which localizes the 
fermions, for one generation, is p2], 

Ydw = hrj^^'^^V + ^5x^5X^5 + /iior,Tr (Wio^io)r? - 2/iioxTr(^ioX^io). (17) 
The resulting 5d Dirac equation, for the charged fermions, is 

--Xiiy)^nYix, y) = 0, n = 5, 10. (18) 

As explained in [22j, to find the zero modes, it is enough to look for solutions for each 
charged fermion of the form \E'„y (x, y) = fnY{y)'ipnY{x) where the ipnvix) are 3+ld massless, 
left-chiral spinor fields. Substituting this into the above Dirac equations yields the solutions 
for the profiles, 

fnviy) = aye-''"-^^), for n = 1, 5, 10, (19) 
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where the Cnv are normahsation constants, and 



Kviy) = hn,j\og (cosh (ky)) + Y\l-hnx arctan f tanh { — 



r _ ^hujv (20) 



h - "'^ 



These profiles have maxima at 



Hence, the charged fermions and the left-chiral neutrino, which reside in the non-trivial 
representations of SU{5), get split along the extra-dimension according to their hypercharges 
and Yukawa couplings to the background. A similar effect was used in [201 El]- 

During the mass fitting sections, we will need to describe the theory in terms of non- 
dimensionalized variables and profiles, since we do not know the value of k. The non- 
dimensionalized domain-wall Yukawa couplings hnx have already been defined in Eq. 
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and so we just need to non-dimensionalise the profiles. Defining the non-dimensionalized 
extra-dimensional coordinate, y, as 

y = ky, (22) 
and changing variables, we see that the normalisation condition for the profiles becomes 

j fny{y)^fny{y) dy = k. (23) 

Hence, in order to use functions which are normalised to one over y, we define the non- 
dimensionalized profiles, fnviy), as 

fnY{y) = k-"^fnY{y). (24) 

Thus, the profiles /„y scale as k^ times a dimensionless function, which is not surprising 
since we know that the 4+ld field \l/„y(a;, y) has mass dimension 2, while the 3+ld field 
'4^nY{x) has mass dimension | as usual. 



Since the factor of k in Eq. 23 will always arise in the normalisation condition when 



changing variables from y to y, we will define the non-dimensionalised profiles for any field 



in the same way as Eq. 24, and they will be denoted with the same symbol used for the 



dimensionful profiles but with an overscript tilde. 



In the case that we have m > 1 generations of fermions, Y^w is generahzed to 



Ydw = hm%v + h'i^nx% + Kl^Trm.^i,),^ - 2hl^Trm,x^{,), (25) 

where i and j are summed from 1 to m. Hence, in the general case, there can be intergenera- 
tional mixing between the quarks and leptons through the interaction with the background. 
The background couphngs and /i„^ have now become 3x3 Hermitian matrices over 
flavour space and need not commute. To solve the equations, we look for zero mode solu- 
tions of the form, 

KYi^,y) = rMyWnYi^), (26) 

where the V'nyl-^) massless left chiral 3+ld fields for n = 5, 10. Putting this into the 
4+ Id Dirac equation results in the matrix differential equation for the profiles /„y, which 
are now 3x3 matrix valued functions of y, 

^^"^^y^^ + v{y)hnvfnY{y) + )J^^xi{y)hnxfnY{y) = 0. (27) 

The case where /i„^ and do not commute, which leads to a natural realisation of the 
twisted split fermion scenario discussed in Refs. [131 E], cannot be solved analytically, and 
so for the sake of simplicity we will only search the parameter space that obeys, 

[hnrj, hnxl = 0, for u = 5, 10. (28) 

Since both the matrices are required to be Hermitian as well, they are thus simultaneously 
diagonalizable, so that for some unitary matrices 5*^, 



(29) 



where the h\^^ and h\^^ are understood to be the eigenvalues of hn-^ and respectively. 
Choosing to localize left-chiral zero modes for is then equivalent to demanding that all 
the eigenvalues of hnr, are positive definite. Solving the 5d Dirac equation then yields the 
general solution for the profiles. 



Uy{y) = S^diag{Clye-''^y^y\ Clye-'ly^y\ ■ ■ ■ , e-''"-(^))Ky, 

(30) 

5-x ^.anhf^l 1 



/3~ f kv 

buYiy) = Kv log (cosh (%)) + Yd -hl^ arctan ( tanh (— ~ 
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Here we have written the multi-generation solutions in terms of the solutions for the one 
generation case. The C^y are normalisation constants, chosen such that the profile matrix 
fnyiy) satisfies the normalisation condition, 

j f\Y{y)fnY{y)dy=l. (31) 

The parameters h}„„ and /i* are the non-dimensionalized versions of and /i* , and are 
defined in the same way as the non-dimensionalized constants from the one generation case 



were in Eq. [20| The V^y are unitary matrices which are present since the solution is unique 
up to matrix multiplication. The VnY-, in fact, correspond to a choice of which 3+ld states 
are the domain wall eigenstates and thus localized to the wall. Unless otherwise stated, we 
will assume these to be the same as the weak interaction eigenstates. 



D. Adding singlet right-handed neutrinos 

To localize the right-chiral neutrinos, we need to couple them to the background. As they 
are gauge singlets, they cannot couple to the adjoint Higgs x- Thus we can only add, 

- h\{WN^r], (32) 

to The relative minus sign in front of the Yukawa interactions for the A^* is introduced 

because for these fields we want localized right-chiral zero modes which represent the right- 
handed neutrinos in the effective 3+ld theory, as opposed to left-chiral zero modes. This 
allows us to treat /ii^ in the same way as /is.^ and /iior?- 

Writing down the 5d Dirac equation for the N\ and demanding that N'^{x,y) = 
f^j^{y)N^{x), where the N^{x) are 3+ld right chiral zero modes, in similar fashion to the 
charged fermions, leads to the profile, 

^1 y ^2 ^ V 

f^(y) = ^m^(C'^e-^ log (cosh (%))^^2g-^ log (cosh (fcy))^... ^ ^nr ^- ^ log (cosh (ky)) ^y^ ^ 

(33) 

where 5*1 is again a choice of basis matrix for the 5d fields, Vn is a change of basis matrix 
for the 4d fields, the C{ are normalisation constants, and the are the positive definite 
eigenvalues of the Yukawa matrix hi^. 
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E. Localizing the electroweak symmetry breaking Higgs boson 



The electroweak breaking Higgs doublet is localized in a very similar manner to the 
fermions. The most general localizing Higgs potential which respects the SU (5) and discrete 
symmetries is 

= + A3($t$)^ + A4$^$?7^ + 2A5$t$Tr[x^] + X^^^ix^^^ + Ay^^x^*?/- (34) 

To find the profiles of the electroweak Higgs doublet, ^^j, and the colored Higgs triplet, $c, 
embedded in the quintet we search for solutions of the form, 

^wA^y y) = PwAy)(pwA^)^ (35) 

where the p^_c are the respective profiles, and 0^^c satisfy the Klein-Gordon equations, 

□ 3+id4>w,c = m^^^0^,c + • • • (36) 

where the masses of the lowest energy modes for ^w,c- Substituting this ansatz into 

the 4+ Id KG equation with the potential V$, one obtains the equations for the profiles 

where 



+ WY{y)PwAy) = ^w.cPwAy)^ (37) 



3y2 /gy 

Wviy) =^^l + W + AsX? + ^^6X? + y 5 Y^7^?Xi- (38) 



Changing variables to the dimensionless coordinate y defined in Eq. [22| the potentials 
of the above Schrodinger equations can be rewritten as shifted hyperbolic Scarf potentials, 
that is we can write them in the form 

+ Al + {Bl- Al- Ay) sech {yf + By{2Ay + 1) sech {y) tanh {y) p^Ay) 

= K,cPwAy) (39) 

where 



-1 + ^2((A5 + ^Ae - A4 - \Y + ^A?)i - 2A5 - ^Ae + 2A4 + ^ 

Ay = 



B 



Y 



(40) 



2((A5 + ^-Ae — A4 — |)2 + %-Af)2 — 2A5 — ^Ae + 2A4 + | 
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and the non-dimensionalized Higgs parameters and masses are defined as 



A4 



fc2 



A ^6 

^6 



(41) 

/^^ — ^2 ' 

2 

~ 2 w,c 

''^■w,c — 

and Xu],c = ^w^c ~ ~ A4 + Ay are the eigenvalues of the equations for the electroweak 
Higgs and the colored Higgs hyperbolic Scarf potentials. 

The hyperbolic Scarf potential has been well studied [28] as it is a member of a class of 
potentials satisfying the shape-invariance condition in supersymmetric quantum mechanics 
(for more on shape- invariant potentials see [221 EO])- For Ay > 0, it is known to have a set 
of discrete bound modes for n = 0, 1, [Ay\ , with eigenvalues 

A:^,, = 2nAy - n\ (42) 

Combining this with the previous equations for A^_c, we see that the potentials have a 
discrete set of bound modes with masses given by 

'^l,u,,c = /i| + A4 - {Ay - nf. (43) 

The physical electroweak Higgs and colored Higgs fields in the effective 4d theory on the 
brane correspond to the n = modes, and they thus exist in the 4d theory if Ay > 0. 
Assuming this, the profiles for these Higgs particles, Pw{y) and Pc{y) respectively, have the 
same form as those of the zero mode profiles for the charged fermions, 

„ _ ri p-bw,c{y) 

bw,c{y) = Ay log (cosh (ky)) + 2By arctan ( tanh ("^) 1 • 

Hence, we can interpret Ay and By to be effective couplings of the Higgs fields to the kink 
and the lump respectively. 
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The effective couplings Ay and By depend on the hypercharges, and thus they are in 
general different for the two Higgs components. This has a number of consequences. Firstly, 
since the masses of the electroweak and colored Higgs depend on their respective Ay, the 
masses of the two components are split. There exists a parameter region where the elec- 
troweak Higgs has a tachyonic mass, < 0, while that for the colored Higgs (if a bound 
state exists) is non-tachyonic, thus inducing electroweak symmetry breaking on the brane 
while preserving SU{3)c, as is desired. Since we know the exact form of the masses, a 
straightforward analysis shows that this parameter region is 

Aly., < fil + \^ < A\. (45) 

Secondly, as there only exist discrete bound modes for a species if Ay > 0, there exist 
parameter regions where the electroweak Higgs component will have discrete bound modes 
localized to the domain wall while at the same time the colored Higgs will have only unbound 
continuum modes in its spectrum. This suggests that an alternate approach to suppressing 
colored Higgs induced proton decay may be possible, as the continuum modes propagate in 
the full 4+ld spacetime so that the partial width contributed to proton decay from these 
modes may be suppressed by further powers of Mqut- The analysis of this situation is 
beyond the scope of this paper. 

Note that it is possible for more than one KK excitation of the Higgs doublet to have 
nonzero vacuum expectation values, thus naturally generating a multi-Higgs doublet model 
on the brane. However, for simplicity, we will choose parameters such that only the elec- 
troweak Higgs has a tachyonic mass, and not its KK modes, and we will also have a bound 
state for the colored Higgs. For the purposes of this paper, we will use three such choices. 

For the first choice, 

/^(j) — k , 
0.5P 



A 



4 



y2 



^5 = ^, (46) 

A7 — — ^) 
vA 

the mass eigenvalues are = — O.SIOA;^, and = 0.380A;^. The graphs of the profiles are 
shown in Fig. |2} 
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FIG. 2: Profiles for tlie colored Higgs, pc, and electroweak Higgs, p^, for the parameter 



choices given in Eq 46 



For the second choice, 



_ 0.5P 

A4 — — 7^ — , 

As 



V 

10A;2 



(47) 



vA ' 

the mass eigenvalues are = — IG.Sfc^, and = 13.2k^. As can be seen in Fig. [sj this 
leads to profiles which are much more localized than those for the first choice of parameters. 
As we will see, this has important consequences for the spread of domain wall Yukawa 
couplings and for the suppression of some of the decay modes for colored Higgs induced 
proton decay. 
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FIG. 3: Higgs profiles for parameter choices given in Eq 
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For the third choice, 



/i| = 97700fc^ 
-75000P 



A4 = 

A5 

A6 = 

A7 



15000P 
-750000P 

2OOOOP 

vA 



(48) 



The resultant squared-masses for the lowest energy modes are = — and = 
2.25 X 10^ A;^. As we can see in the graphs of the profiles in Fig. El for this parameter region, 



3.0 




2.5 




2.0 




1.5 




:' 




;' '. 0.5 





P,(.v) 



FIG. 4: Higgs profiles for parameter choices given in Eq 
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the electroweak Higgs is highly peaked near the brane at ?/ = 0, while the colored Higgs is 
more delocalized and substantially displaced from the wall. This parameter choice exploits 
the property of the Higgs sector that effective kink and lump couplings Ay and By are 
not the same for the colored and electroweak Higgs. As we will see, this kind of parameter 
choice can lead to suppression of all decay modes for colored Higgs induced proton decay, 
and ensure that the partial lifetimes for these modes are all many orders of magnitude above 
the current lower bounds. 

Note that the Higgs vacuum expectation value, {4>w), is not uniquely determined by the 
constants which determine the Higgs profile. By dimensional reduction of the action, 

S = j d^x - 14, (49) 

one can show that the effective electroweak symmetry breaking potential is 

Vew{(I)w) = A'(0]^0^)^ + ml(j)l(f)y„ (50) 

16 



where 

A' = As I ptiy) dy. (51) 
Thus the VEV of the Higgs doublet is 



—rn^.. 



-ml 



2X' 

(52) 



2\sfpi{y)dy' 

and so whatever we choose for the other constants, we can always adjust As appropriately 
so that we get the correct VEV of 174 GeV. 



F. Generating mass matrices for the charged fermions 

The electroweak Yukawa Lagrangian, Y5, from [22] which generates masses for the charged 
fermions is generalized to 



Y, = h^{%)c^i^{,^ + hXe-^^'-i<i^\,)0^,^,^sK + h.c, (53) 

for m generations of fermions. Here, lower case Greek letters are SU{5) indices and the 
lower case Roman letters indicate flavor. 

The /i_ terms generate mass matrices for the down-type quarks and electron type leptons, 
while the h+ terms generate a mass matrix for the up-type quarks. Extracting the compo- 
nents from each term which generate 3-|-ld masses and performing dimensional reduction, 
one finds the mass matrices to be 

M^ = Av flJy)h+fQ{y)p.M dy, 



Ma = ^v j fl^{y)h_fQ{y)pM dy, (54) 

Me = -^v j flj^{y)h^fi{y)Pw{y) dy, 

where Pwiy) is the profile of the electroweak Higgs doublet which is embedded in $, and 
V = 174 GeV is the vacuum expectation value of the electroweak Higgs field attained on the 
brane. 

Converting to dimensionless quantities, and defining the non-dimensionalized electroweak 
Yukawa couplings by 

h+ = kh+, 

(55) 
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we see that these mass matrices can be rewritten as 




j flj,iy)h+fQ{y)Pw{y) dy, 
^ j fdR(y)^-fQ{y)PM dy, 
^ / fLiy)h-fiiy)Pu,{y) dy, 



(56) 



There are some important consequences of the above forms of the mass matrices, which 
depend on overlap integrals of the profiles for the left and right chiral fermions and the elec- 
troweak Higgs. Firstly, the overlap integral dependence means we avoid the usual incorrect 
mass relations like rrie = rrid which are characteristic of ordinary 3+ld SU{5) models with 
a Higgs quintet. This is also the reason why we do not need a Higgs belonging to the 45 
representation of SU (5) containing an electroweak Higgs triplet to get the Georgi-Jarlskog 
relations |3T]. Thirdly, since the fermions are split according to their hypercharges, and 
the splittings are dependent on the background couplings, we can potentially generate the 
fermion mass hierarchy and mixings by splitting the fermions appropriately so that the over- 
lap integrals are in the desired ratios. It will be shown in a later section that this can be 
done. 

G. Generating Dirac neutrino masses 

To generate Dirac masses for the neutrinos, we need to add Yukawa interactions involving 
the "$1, which contain the left handed neutrinos, the A^*, which contain the right handed 
neutrinos, and $ which contains the electroweak Higgs. The correct terms to add to Y5 
which are both SU{5) invariant and respect the reflection symmetry which preserves the 
topological stability of the domain wall are 



Reducing these terms to their SM components, and integrating out the extra-dimensional 
dependence, one finds the resulting Dirac mass matrix for the neutrinos to be 



{hly^%<i>N^ + h.c. 



(57) 




(58) 



Defining the dimensionless neutrino Yukawa couplings as 



/is = k^h 



(59) 
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and changing to non-dimensionalized quantities, we can rewrite the Dirac mass matrix for 
the neutrino as 

m^^v J fl^{y)hji{y)p^{y) dy. (60) 
H. Generating Majorana neutrino masses 

Let us consider one generation first. To generate a Majorana mass for the neutrino, we 
need to add terms to the Lagrangian that will dimensionally reduce to terms proportional 
to i^^i^i? in the effective 4d theory. Thus, we might want to consider adding a term like 

NN^ + h.c. (61) 

This is obviously gauge invariant, and it turns out that it is also invariant under the discrete 

reflection symmetry as well. Wc first need to consider what implications the addition of this 
term has for the existence of solutions of the 5d Dirac equation. The relevant Lagrangian is 

Ln,dw = iNV^dMN + hx^NN-q - \m{NN^ + WN), (62) 

and thus the 5d Dirac equation becomes 

iV^dMN + hir,N7] - mN^ = 0. (63) 

Demanding the conditions that 

N{x,y) = fN{y)'^R{x), 

l^T^R = (64) 

and noting that the parts proportional to and {i'rY must be independent of each other 
as the corresponding spinors transform as right-chiral and left-chiral spinors respectively, we 
get two independent equations for /at, 

^ + hirjri{y)fN{y) = 0, 

(65) 

m'fNiy) - imf^{y) = 0. 
The first of the equations above is exactly the same differential equation as before without 
the new term, and thus the /jv must also have the same form as before, 

f^{y) = C-Are-"^ '°s(cosh(fcj/)). (66) 
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The second condition then imphes that m' = \m\, and since any phase can just be absorbed 
into the definition of A^, we can take m' = m. Hence, instead of a right-chiral zero mode, 
we now have a right-chiral Majorana mode of mass m locahzed to the domain wall. 

Similarly with three generations, the profiles are unaltered by the Majorana mass terms, 
and the 3+ld Majorana mass matrix after dimensional reduction is then 



We have thus successfully shown that both Dirac and Majorana masses can be generated 
with the addition of a right chiral singlet neutrino, and thus the see-saw mechanism can 
be employed. We will now demonstrate that the fermion mass hierarchy and small CKM 
mixing angles can be generated from split fermion idea. 

III. GENERATING THE FLAVOR HIERARCHY AND MIXING ANGLES 

The fermion mass matrices depend on overlap integrals of the fermion profiles and the 
electroweak Higgs. Since the left-chiral and right-chiral components are naturally split ac- 
cording to their hypercharges, and since these overlap integrals are exponentially sensitive to 
these splittings, it seems we can employ the split fermion idea to account for the fermion 
mass hierarchy from a set of domain wall couplings which are all about the same order of 
magnitude in this model. 

Throughout the rest of the paper, we will quote the dimensionless background Yukawa 
couplings to five significant figures. The reasons for this are the exponential sensitivity of 
the profiles to these couplings and the difficulty that was found in generating the neutrino 
mass squared differences (which are quadratic in overlap integrals of these profiles) to an 
acceptable and reasonable precision. Since this is also a classical calculation where quan- 
tum corrections are ignored, and since the quark and neutrino masses are not as precisely 
measured or well known as those for the charged leptons, we will quote the resultant masses 
of the quarks and neutrinos to two significant figures, neutrino mass squared differences to 
one significant figure, and the charged lepton masses to three significant figures. 




(67) 
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A. The one-generation case with a Dirac neutrino and the suppression of colored- 
Higgs-induced proton decay 



In this section we shall show that the mass hierarchy amongst the first generation of 
fermions can be generated from the split fermion idea [9] which arises naturally in our 



model. We will start with looking for solutions with the Higgs parameter choices of Eq. 46 



Firstly, we must make the neutrino light. The right chiral neutrino is always localized 



at y = while the choice of Higgs parameters in Eq. 46 (and in fact for those in Eqs. 47 



and 48 as well), the Higgs is localized to the right. Hence, the easiest way to induce a small 
Dirac neutrino mass is to shift the lepton doublet to the left. As the lepton doublet, L, 
has hypercharge —1 and the charge conjugate of dn has hypercharge +|, choosing ^ to 
be negative will displace the lepton doublet as desired while placing the right-chiral down 
quark to the right, near the electroweak Higgs. We now need to make the charged fermion 




FIG. 5: The profiles for ur, L, (Ir and the electroweak Higgs with the Higgs parameter 
choice of Eq. 46, /ii^ 



115, hr, = 100 and k 



5x 



-250. 



masses significantly larger. Since the charge conjugates of ur and en, and the quark doublet 
have hypercharges — |, +2 and +^ respectively, making the ratio positive will shift 
far to the left, towards the lepton doublet, Q to slightly to the left, and ur to the right. 
We found the following solution by using this configuration, making the parameter choice 
h^rf = 100 and h^^ = —250, plotting the contours along which the overlap integrals give the 
desired mass ratios, and then finding where the two contours intersected. Doing this yielded 
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FIG. 6: The profiles for the 10 representation, and the electroweak Higgs with the Higgs 
parameter choice of Eq. 
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hwrj = 8.2674 and h^x = 27.911. 




FIG. 7: The profiles for all fermions and (j)w for the first solution with hi„ = 115, h^n = 100, 



h 



5x 



250, hio„ = 8.2674, hioy = 27.911 and the Higgs parameter choice of Eq. 



46 



the solution for the couplings for the 10 multiplet, /iior? = 8.2674 and hio^ = 27.911. With 
the ratios now fixed, setting the 5d electroweak Yukawas = hj^ = = 5.2268 x 10^^k~^ , 
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and setting the kink coupling for the right handed neutrino to hi^j = 115 gives the masses, 

= 0.13 eV 
me = 0.511 MeV 

(68) 

fUu = 2.5 MeV 
rud = 5.0 MeV. 

Thus, we have generated an neutrino mass below the current most stringent upper bounds 
of roughly 2 eV [32] , the correct electron mass, and up and down quark masses within current 
constraints of 1.5 MeV< m„ < 3.3 MeV, and 3.5 MeV< < 6.0 MeV [32]. 

Furthermore, it turns out we get significant suppression of some modes of colored Higgs- 
induced proton decay with this setup. The colored Higgs scalar can induce the decays 
p —7- e+vr'^ and p — t- I'e'^'^, for which the Feynman diagrams are shown in Figures [S^a) and 
^h) respectively. 
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FIG. 8: Feynman diagrams for the processes (a) p — )■ e'^n^ and (b) p — )■ I'eTr'^ 
For the process p e'^n'^, the partial lifetime of each contribution is 

(69) 



4 
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where C^^^ and C„^^ are replaced by the effective 4d couphngs strengths of the operators 
inducing the vertices Vue and Vud respectively. The operators responsible for the vertex Vue 



are {eR)'^UR(j)c. and LQ(j)c, and their respective coupling strengths are 

'^J^^un^. = I fnRi.y)fen{y)Pc{.y) dy, 



^LQ^c " 7f ^- / fL{y)fQ{y)pc{y) dy. 



(70) 



The operators responsible for the vertex Vud are {ujiYdRCpl and e^-'''QiQj{4>l)k, and the asso- 
ciated coupling strengths are 



C?^d„<6t = \,h- I fuAy)fdn{y)pc{y) dy, 



{uB.YdR<i>i 

\2, 



(71) 



Cqq,^* = 4/i+ / {fQ{y)yVc{y) dy. 



Similarly, the partial lifetime of each contribution to p — t- z/evr"*" is 

4 



(72) 



The operators responsible for the vertex Vud in the p — t- z/gTr^ are the same as that for 



p — )■ e'^vr", while the operators responsible for the Vdu vertex are LQ^c and {i'RYdB,(t>c- The 



coupling strength for the operator {i'Rydji(j)c is 

^I^dR4,, = ^3 j fun{y)fdniy)Pciy) dy. (73) 

For the solution given above, it turns out that the partial width for p — )■ e+7r° involving 
just right chiral fermions is substantially suppressed, with Cj^^^^^^ = 9.6 x 10~^^ and 
^{unY dn^* ~ '^■^ ^ 10^^. Since the partial lifetime for p — )■ e^Ti^ is at least 8.2 x 10^^ years |33] . 
and given nip = 0.938 GeV [32J, this sets a lower bound on the colored Higgs mass of about 
3.3 X 10^ TeV, much reduced compared to the standard result of rric ~ Aqut ~ 10^^ GeV. 
Since Cj^j^^^^^ = 3.7 x 10~^, and the lower bound of the partial lifetime for p — )■ un^ is 
2.5 X 10^^ years [32|, the contribution to p — ?■ Ueir^ involving the vertices {i'RYdR(f)c and 
(uR)'^dR(f)* sets a lower bound on the colored Higgs mass of 4.8 x 10^ TeV. 

However the contribution involving just left-chiral fermions is not substantially suppressed 
from the splittings, with Cj^q^^ = 9.1 x 10""^ and Cqq^* = 1.0 x 10^^. These operators 
contribute to both p — )■ e~*"7r° and p — )■ z/eTr"*", and thus the partial widths coming from the 
combination of these operators set lower limits on the colored Higgs mass of 7.0 x 10^^ GeV 
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for p — > e~^7r^, and 1.6 x 10^^ GeV to suppress p — )■ z/evr"*". Therefore, we must still fine-tune 
so that the colored Higgs mass is of the order ~ 10^"^ GeV to suppress all proton decay 
modes induced by the colored Higgs scalar. 

One might then ask how to suppress proton decay even further. We could try looking for 
a solution where the profiles are more spread out in the extra dimension. It turns out the 
choice of parameters, hi^j = 100, h^j] = 100, h^^ = —700, /iior? = 0.81688, hio^ = 23.868, 
and = h_ = h-^ = 0.11177A;~^, yields the same masses for the electron and the quarks as 
the first solution, and gives a neutrino mass of the order 

~ 10-2^ eV. (74) 
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FIG. 9: The spread of profiles for the second solution with hi„ = 100, h. 
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100, 



= -700, hior, 



0.81688, /iiox 



23.868, and the Higgs parameter choice of Eq. 



46 



For proton decay, we now have for the operators involving just right chiral fermions 
C„ = 4.2 X 10-l^ Ct-^, = 1.7 X 10-^ and Cj-^, , = 2.8 x lO^^s which 
yield nic > 6.1 x 10^ TeV from the partial width for p — )• e^vr'', and nic > 1.2 x 10^ GeV 
from the partial width contributed to p — > Ueir^. However the couplings for the operators 
involving the left chiral fermions are still not suppressed enough to solve the doublet-triplet 
splitting problem naturally, with Cj^q^^ = 2.0 x 10^^ and Cqq^* = 0.21, setting the bound 
rric > 1.5 X 10^^ GeV from the more constraining decay p — t- e^n^. 
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By spreading out the profiles, we have increased the spread of the domain wall parameters 
while only suppressing the proton decay modes induced from right chiral fermions by a 
further two orders of magnitude. It turns out that a choice of Higgs potential parameters 
giving more localized Higgs profiles can solve the first problem while yielding a similar result 
for proton decay. A solution for the second Higgs profile for which the Higgs parameters are 
is h+ = = = 82.975A;-i, hi^ = 200, h„ = 100, = -250, 



those in Eq. 
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*5x 



60.126, hiQ^ = 99.829, which again yields the same masses for the electron and the quarks 
as the previous solutions, and the neutrino mass 



m,. 



0.024 eV. 



(75) 



Interestingly, for these parameters, C-. 



2.0 X 10-^4 and C-. 



{uRYdiicl)l 



1.1 X 10"^ 



suppressing the mode of p — t- e+7r° involving just the right chiral fermions to the extent 
that the lower bound for the colored Higgs mass set by this mode is just 1.1 TeV. However, 
the decays involving the left chiral fermions are in fact enhanced rather than suppressed 



by the fermion splittings with C 



2.0 X 10-2 and C, 



QQ4''c 



39, which means that the 



suppression factor coming from the effective coupling constants is of order 1, and so we 
need to make ttIc ~ Mqut- Also, the partial width involving the right chiral fermions for 
p — > z/e7r+ is not as suppressed this time, with Cjjj^^^^^ = 8.5 x 10"'', so that an of order 
10^^ GeV is required to suppress this particular mode. 




FIG. 10: The profiles for the solution with hi„ = 200, h^r. = 100, h^^ = -250, 



/iior? = 60.126, /iiox 



99.829, and the choice of Higgs parameters in Eq. 
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The ultimate reason we have successfully suppressed the modes of proton decay involving 
just the right chiral fermions but not those involving the left chiral fermions so far was that 
the vertices involving the right chiral fermions depended on the profiles for ur and dji which 
were localized near the electroweak Higgs away from the colored Higgs, whereas due to the 
setup to generate the mass hierarchy, the quark and lepton doublets were placed significantly 
closer to the colored Higgs. To keep the natural solution to the mass hierarchy problem, we 
do not wish to displace the quark and lepton doublets; a more fruitful option is to choose 
Higgs parameters such that the colored Higgs is well displaced from the domain wall, while 



at the same time the electroweak Higgs is close to y = 0. We have seen in Sec. HE that 
this is in fact possible with the Higgs parameter choice given in Eq. |48j A solution for this 
third Higgs profile to the mass hierarchy problem is the parameter choice h+ = = = 
40987A;-i, hir, = 1000, hr, = 1000, = -1000, /iio^ = 624.62, hio^ = 382.43, which 
yields the same electron, up and down quark masses as before and a neutrino mass of the 
order 

~ 10"^ eV. (76) 

-139 



This time, for the proton decay inducing interactions, Cj—^ , ~ 10 , C- 



^ ^ , 1 C\ 126 ^r\-v 4-TiCi /~\T-\CiT*n -i~/~\T*ci i m r/~\l t nn c 

-92 



10 and C-, — Yd ni ~ ^0 f*^^ operators involving just right chiral fermions and 



for those involving the left chiral fermions, Cj^q^^ ~ 10^ and Cqq^* ~ 10~ . Hence, 
all the decay modes are suppressed by roughly 90-100 orders of magnitude, with the most 
constraining decay mode p — t- e+vr^ with the left chiral fermions now setting a lower bound 
on the colored Higgs mass of ~ 10~^^ eV. Realistically, for such a solution, the colored Higgs 
mass should still at the very least be 45 GeV since we have not seen the Z boson decay into 
them, and more probably ~ 1 TeV since it is proportional to k in this model, so the partial 
lifetime arising from colored Higgs induced proton decay would be over 10^^^ years. 

The profiles for the latter two solutions are much less dispersed than those first two for 
the most delocalized electroweak Higgs profile, as one might expect. This is reflected in the 
breadth of the domain wall parameters; for the third and fourth solutions, the ratios of the 
magnitude of the largest background coupling {h^^ = —250 and /ii^ = h^rj = —h^x ~ 1000 
respectively) to the smallest (/iio»; = 60.126 and /iio^ = 382.43) are roughly 4.2 and 2.6 
respectively, under an order of magnitude. In comparison, for the first two solutions, the 
corresponding ratios are about 30.2 and 860 respectively. This difference, as we will see, is 
exacerbated for three generations. 
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FIG. 11: The profiles for tlie solution with hn = 1000, hr, = 1000,/i5x = -1000, 
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624.62, h 



382.43, and the choice of Higgs parameters in Eq. 
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In summary, we have shown that the one generation mass hierarchy can be generated by 
splitting the fermions without fine tuning the electroweak Yukawa constants, and that for 
an appropriate choice of Higgs parameters, one can reduce the spread of the domain wall 
Yukawa constants and suppress proton decay by roughly 100 orders of magnitude without 
fine tuning the colored Higgs mass. In the next section, we will see that we can also do this 
for the three generation case without quark and lepton mixing. 

B. Generating the higher generation mass hierarchies without electroweak mixing 

For the sake of simplicity, we can solve the mass hierarchy problem while first omitting 
quark and lepton mixing by setting the off diagonal elements of the electroweak Yukawa 
matrices to zero. Solutions are found in analogous fashion to the one generation case by 
finding where the overlap integrals are in the desired ratios in domain wall parameter space. 

For the first Higgs profile, setting h'^_^ = = h\ = lA093k~^ , the set of parameters 
choices for the domain wall Yukawa parameters and the resultant masses for the electron- 
type leptons, m\j for the up-type quarks and m}, for the down-type quarks are shown in 
Table |ll As one can see, these masses all lie within current experimental limits [32]. We 
then get similar results for the second Higgs profile with h\ = = h\ = 7859. 3A;~^ and 
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the third Higgs profile with hi 



hi 



2 70 1.2 A; 2 with the solutions for these two 



parameter choices given in Tables M and III respectively. 



















i 


'''bri 




'nori 




m'^(MeV) 


mi^(MeV) 


m*o(MeV) 


1 


1064.0 


-8563.9 


0.2 


25.496 


0.511 


2.5 


5.0 


2 


48.986 


-708.28 


1.5 


17.330 


106 


1.3x10^ 


1.0x10^ 


3 


100 


-300 


10.537 


8.6032 


1.78x10^ 


1.7x10^ 


4.2x10^ 



TABLE I: A set of domain wall parameters and the resultant masses with Higgs parameters 
chosen in Eq. 
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and electroweak Yukawas set to h'. 



hi = h\ = 1.410A;-2 for i = 1, 2, 3 



















i 








'T'lOx 






m\j{MeY) 


1 


200 


-648.41 


38.552 


99.220 


0.511 


2.5 


5.0 


2 


200 


-493.42 


62.128 


94.251 


106 


1.3x10^ 


1.0x10^ 


3 


200 


-400 


73.744 


76.383 


1.78x10^ 


1.7x10^ 


4.2x10^ 



TABLE II: A set of domain wall parameters and the resultant masses with Higgs 
parameters chosen in Eq. 47 and electroweak Yukawas set to hi 

for i 



W = hi 



7859.3fc-2 



1,2,3 



There exists a finite range of parameter space spanned by the remaining couplings h\^ 
of the right handed neutrinos to the domain wall which fit the currently accepted squared 
neutrino mass differences of Amu = T.Q^gl x 10"^ eV^ El and Am23 = 2.74;^-26 x lO'^ eV^ 



[35] . and cosmological constraints (in some models) on the sum of the masses < 
0.3 — 0.6 eV [36| 137] , for normal, inverted, and quasidegenerate neutrino mass hierarchies. 
Provided that the fermion doublets are sufficiently localized and displaced away from the 
domain wall at y = 0, where the right handed neutrinos are always situated, and the 
electroweak Higgs, one can just adjust the couplings of the right handed neutrinos to the kink 
to get the desired masses and hierarchy. For each of the three solutions given in Tables [1} [TT] 
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i 


"'5n 


""Sy 


'^10?? 


"■lOv 


m*^(MeV) 


mi;(MeV) 


m*^(MeV) 


1 


2000 


-1585.2 


660.91 


369.07 


0.511 


2.5 


5.0 


2 


2000 


-1434.5 


744.05 


325.26 


106 


1.3x10^ 


1.1x10^ 


3 


2000 


-1300 


708.14 


256.02 


1.78x10^ 


1.7x10^ 


4.2x10^ 



TABLE III: A set of domain wall parameters and the resultant masses with Higgs 

and electroweak Yukawas set to h\_ = W_ = h\ = 2701.2A;~5 
for i = 1, 2, 3 



parameters chosen in Eq. 



48 



and 



III 



three example parameter choices for the h\ yielding normal(N), quasidegenerate(Q), 



and inverted(I) neutrino mass hierarchies are given in Tables IV, Vjand VI respectively. 



Hierarchy 






"It? 




m^^(eV) 


m^^ (eV) 


N 


100 


15.44 


110.4 


1.5x10-"^^ 


0.0089 


0.053 


Q 


18.919 


13.764 


106.61 


0.10 


0.10 


0.086 


I 


19.503 


14.219 


300 


0.051 


0.052 


7.4x10"^ 



X;?Tii(eV) 


Am2^(eV2) 




0.062 


8x10"^ 


3x10-3 


0.29 


8x10-5 


-3x10-3 


0.10 


8x10-5 


-3x10-3 



TABLE IV: Solutions for normal, quasidegenerate and inverted neutrino mass hierarchies 

given the parameter choices given in Table jlj 

The distribution of the fermions for each family for both the solutions with a normal 



neutrino mass hierarchy are shown in the Figs. 12, 13 and 14 As can be seen, the lighter 
generations are, on average, more spread apart, more distant from y = and more delocal- 



ized. Comparing the plots for the first Higgs profile in Fig. 12 to those for the second and 



third Higgs profiles in Figs. [13] and [T4| it is conspicuous that this increase in spread of the 
fermions between generations is dramatically reduced for the more localized Higgs. This is 
reflected in the spread of domain wall parameters, with the ratios between the smallest and 
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Hierarchy 






"■If] 


m^^(eV) 


m^^(eV) 


rriu^ (eV) 


N 


200 


132.73 


262.60 


3.5x10-15 


0.0089 


0.053 


Q 


54.564 


114.28 


253.20 


0.096 


0.096 


0.081 


I 


56.690 


118.95 


650 


0.051 


0.052 


2.9x10-6 





Amii(eV2) 




0.062 


8x10^5 


3x10-3 


0.27 


8x10-5 


-3x10-3 


0.10 


8x10-5 


-3x10-3 



TABLE V: Solutions for normal, quasidegenerate and inverted neutrino mass hierarchies 

given the parameter choices given in Table |TT| 



Hierarchy 


hi 


'hri 






m^^(eV) 


niv^ (eV) 




Amii(eV2) 


Ami2(eV2) 


N 


2000 


1449.2 


2044.3 


8.1x10-^ 


0.0089 


0.053 


0.062 


8x10-5 


3x10-3 


Q 


826.28 


1250 


1948.5 


0.084 


0.084 


0.099 


0.27 


8x10-5 


3x10-3 


I 


852.44 


1291 


4500 


0.051 


0.052 


2.3x10-6 


0.10 


8x10-5 


-3x10-3 



TABLE VI: Solutions for normal, quasidegenerate and inverted neutrino mass hierarchies 

given the parameter choices given in Table 



III 



largest non-dimensionalized domain wall parameters. For the parameter choices of Tables |T], 



III and III and normal neutrino mass heirarchies, these ratios are respectively 4.3 x 10 , 17, 



and 7.8. The solution of Table III is particularly interesting since the non-dimensionalized 
electroweak Yukawa constant /?■!__,_ 3 = h'^_^^k^ = 2701.2 is of the same order as the non- 
dimensionalized domain wall parameters. 
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FIG. 12: Plots of the profiles of the first (a), second (b), and third generation (c) of fermions 
with the parameter choice of Table [T] and the of the normal hierarchy parameter choice in 
Table [IV 
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FIG. 13: Plots of the profiles of the first (a), second (b), and third generation (c) of fermions 
with the parameter choice of Table [TT] and the of the normal neutrino mass hierarchy pa- 
rameter choice in Table fVl 
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FIG. 14: Plots of the profiles of the first (a), second (b), and third generation (c) of fermions 



with the parameter choice of Table |III| and the normal hierarchy parameter choice in Table 

vn 
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With regards to proton decay, the results for the parameter choices of Tables |H |TT] and [ITT 



are similar to those of Sec. Ill A[ For the first two solutions with Higgs parameters chosen 



from Eqs. 46 and 47, the decay modes involving just right-chiral fermions are substan- 
tially suppressed while there is negligible suppression for the modes involving the left-chiral 
fermions. For the parameters chosen in Table [l| Cj^r^^^^^ ~ 10~^^, Cj^^^^^^, = 1.7 x 10~^, 
^LQ<f>c ~ ^•'^ ^ ^^^^ ^QQ<t>i ~ 0.80, and hence the lower bound on the colored Higgs mass 
from the decay mode for p — t- e+7r° involving just right-chiral fermions is roughly 120 TeV, 
while that from the decay mode involving just left-chiral fermions is of order 10^^ GeV. 



For the solution in Table 



II 



, ~ lO-'^", , = 4.1 X 10-^ CjoA, = 2.5 and 

CqQff,* = 7.4 X 10^, so that the decay mode involving the right-chiral fermions sets a lower 
bound on rric of just 1.7 GeV, while the decay mode involving just left-chiral fermions are 
in fact enhanced, with the lower bound on rrif. increased to order 10^'' GeV. 



For the solution of Table III, just as with the one generation solution with the Higgs 



parameter choices of Eq. 48, all decay modes are suppressed since the colored Higgs is 
well away from the domain wall and the electroweak Higgs. For this solution, we have 
Cj^u,,. - 10-l3^ Cj^,^,, ~ 10-l2^ Cj:q,^ ~ 10-9^ and Cqq^, ~ IQ-^', so that the 
bound on rric set by the less suppressed decay mode involving the left-chiral fermions is of 
the order of 10"''^ eV. For all neutrino mass hierachies, the coupling constant C-^jj^^^^ is 
also well below 10^^°'^, so that p — v^.t^'^ is also negligible. 

Now that it has been demonstrated that the three generation mass hierarchies can be 
generated from the exponential dependences of the overlaps on the domain wall couplings 
while suppressing proton decay, the next step is to incorporate quark and lepton mixing. 



C. Accounting for the Cabibbo angle in the two-generation case 

To produce realistic mass matrices, we must account for the fermion mixing angles as 
well as the masses. For the sake of simplicity, we will work with two generations and show 
that the Cabibbo mixing angle can be produced along with the mass hierarchy. 

Performing the required analysis is quite complicated since, if we are to assume that 
all the 5d Yukawa couplings are equal, including the off diagonal couplings, then the order 
of the analogous equations giving desired mass matrix element ratios from the overlaps is 
the equal to the number of families in the theory. Thus, to generate the Cabibbo angle, 



35 



one must solve equations which are quadratic in the overlaps, and which are also no longer 
separated with respect to the domain wall parameters. For the CKM matrix it is even 
worse since the equations are cubic. This raises difficulties, in particular, with the down 
and electron sectors, since these sectors depend on all of the background couplings of the 
charged fermions, which amount to eight for two generations. Hence, we are forced to start 
with the up quark sector ffist, for convenience, since it only depends on four couplings. This 
makes it difficult to guarantee that the Dirac neutrino masses will be light. 

Instead of directly solving the equations quadratic in the overlaps, we will try to generate 
mass matrices approximately equal to the Cholesky decompositions of the desired mass 
matrices squared, M^M. This is similar to an approach of generating mass matrices in 
NNI (Nearest-Neighbour-Interaction) basis for the three generation case, as was done in 
the analysis with Gaussian profiles done in |38], and by Mirabelli and Schmaltz [39]. The 
advantage of this approach is that we can now do the analysis in terms of equations linear in 
the overlaps instead. The main disadvantage is that we must rely on one of the off diagonal 
terms being significantly suppressed compared to the other couplings. 

It turns out that this approach can get the charged fermion mass hierarchies and the 
Cabbibo angle. Taking all the electroweak Yukawa couplings to be = = /ig = 

and making the choices for the 



0.089104/c 2, the Higgs parameters to be those of Eq. 
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domain wall Yukawas in Table VII, we obtain the following mass matrices correct to three 













i 


'''5ri 




'nori 




1 


12.585 


-36.719 


100 


53.346 


2 


365.78 


-1708.2 


2.5273 


27.095 



TABLE VII: The choices for the domain wall parameters with Higgs parameters from Eq. 

and = h'[ = lii = 0.089104fc-i for i,j = 1,2,3 
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significant figures. 



1.04 X 10^ 7.28 X 10^ 
2.77 X 10-3 3.05 



(77) 
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70.3 77.8 
'0 7.47, 



(78) 



and 



(79) 



106 ~ 
0.465 0.511 

The resultant left diagonalization angle for M„ is 0^ = 55°, and the left diagonalization 
angle for Ma is 9^ = 42°, yielding the Cabibbo angle 9c = — = 13°. Taking the 
square roots of the eigenvalues of M^M„, M^M^, and MjMg yields the masses of the charged 



fermions, and they turn out to be those in Table VIII 



i 


m^(eV) 


mlj{eV) 


m})(eV) 


1 


0.511 


2.5 


5.0 


2 


106 


1.3x10^ 


1.1x10^ 



TABLE VIII: The masses of the mass eigenstates in the electron, up and down type 
sectors respectively with the parameter choice of Table 



VII 



This solution does not permit two light neutrino masses for the case of a Dirac neutrino. 
For example, for the parameter choice h\ = h\ = 100, the mass matrix for the neutrino is 



6.23 ~ 
6.23 ~ 



(80) 



which yields mj,i = and 171^2 = 8.8 MeV. Because for this solution one of the left weak 
eigenstates of the neutrino has a h\^/h\^ ratio larger in magnitude, and thus is more distant 
from the wall and the Higgs, the mass matrix for the neutrino for this solution will always 
have the entries of one column being larger than the other for significantly large and 
since one of the eigenvalues is typically of the same order as the larger of the two elements 
in the larger column in such a matrix. In the limit — )■ 00 for i = 1, 2, 



(81) 



M'J = hv J U^{y)fL^{y)pM dy 
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as the profiles for tlie riglit lianded neutrinos converge to delta functions at y = 0, one 
obtains the mass matrix, 

/5.81 ~o\ 

M,= \ , (82) 

\^5.8i ~oy 

which yields neutrino masses of mi = 8.2 MeV and 1712 = 0. Thus, the best we can do for 
this solution which yields Cabibbo mixing and the charged mass hierarchy, for the case of 
a Dirac neutrino, is to generate a massless neutrino, and a neutrino about 2-3 MeV heavier 
than a down quark. 

This does not prove that there is no solution for a Dirac neutrino which incorporates 
quark mixing and the fermion mass hierarchy. The scheme we used and the section of 
parameter space searched led to one of the lepton doublets being placed too close to the 
right handed neutrino and too delocalized to support two light neutrinos. A more thorough 
search of the parameter space, perhaps utilizing a Monte Carlo method, will have to be 
done to determine whether a solution supporting two sufficiently light neutrinos exists for 
the case of a Dirac neutrino. 

For a Majorana neutrino, however, this solution presents no such problems. As was shown 
earlier in this paper, the seesaw mechanism can be employed in the model, and can thus be 
used to suppress the mass of the heavier neutrino. 



The set of domain wall parameters in Table |VII| generate the desired mass spectrum and 
the Cabibbo angle. The ratio between the parameters smallest {HIq^ = 2.5273) and largest 
{hl^ = —1708.2) in magnitude is roughly 670. To reduce this, we would need to find solutions 
with a more localized electroweak Higgs, such as those resulting from the parameter choices 



Eqs. 47 and 48 Finding such solutions has been difficult, however, and we will leave this 
to be done in a Monte Carlo search. 

Proton decay for this solution cannot be suppressed without fine-tuning the colored Higgs 



mass. Since we are using the Higgs parameter choice of Eq. and not that of Eq. 48 



the colored Higgs is sufficiently close to the domain wall so that the decay modes involving 
just the left-chiral fermions are not sufficiently suppressed. Furthermore, since the off- 
diagonal electroweak Yukawa constants are now non-zero, operators such Uji{^f()(f)* and 
~SR{uRY(j)c are present in the action. This means we also have to account for the decay 
modes p — ?■ /i~^7r°, p — )■ e~^K^ and p — ?■ ^'^K^, which have partial lifetime lower bounds of 
6.6 X 10^^ years [33j, 1.5 x 10^^ years [32j and 1.2 x 10^^ years [32j. After computing the 
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overlaps to find the interaction strenghts in the weak eigenbasis, transforming to the mass 
eigenbasis and then performing similar analyses for each of the decay modes, we find that 
even the most constraining decay mode involving just right-chiral fermions, that of — > /x^vr^ 
sets a lower bound on the colored Higgs mass of 3.2 x 10^ GeV. In fact both the modes 
involving the antimuon set higher bounds then those producing positrons since the coupling 
for the ui^{fi[i)(j)* vertex is a few orders of magnitude higher than that for the UR{eii)(f)* 
vertex, due to right handed muon being closer on average to the wall. The decay modes 
involving left-chiral fermions are still barely suppressed, with the decay modes involving just 
left-chiral fermions for both p — )■ e+7r° and p — ?■ fi^ir^ setting lower bounds on rric of order 
10^^ GeV. Obviously, this situation would change if we used the Higgs parameter choice of 



Eq. 48 for which the colored Higgs is well displaced from all fermions and the domain wall. 



D. Lepton mixing 

It appears that generating near tribimaximal mixing in the lepton sector is incompatible 
with the results for the fermion mass hierarchy problem, for both Dirac and Majorana 
neutrinos. 

As we have seen, solutions to the mass hierarchy problem typically involve shifting the 
lepton doublets to different locations away from the domain wall and the electroweak Higgs. 
This means that for such solutions, assuming the domain wall couplings for the right handed 
neutrinos are roughly equal, the neutrino mass matrix will take the form. 



-2 



(83) 



where e ^ 1 and the constants Oj, bi and Cj are all taken to be roughly the same order of 
magnitude. Then the mass matrix squared will take the form 



^ |ap a. be a.ce^^ 
a.be |bpe^ b.ce^ 
a.ce^ b.ce^ IcPe'^ 



(84) 



MIMu is clearly hierarchical, and thus the neutrino sector cannot generate two large 
mixing angles. From the electron sector, we know that both the lepton doublets and the 
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right handed electrons are placed away from the electroweak Higgs doublet, and thus their 
overlaps decrease rapidly with the splittings, inducing hierarchical electron mass matrices. 
Hence, we cannot generate large mixing angles in the electron sector either, with the generic 
type of solution for the mass hierarchy, and thus tribimaximal mixing cannot be produced 
in the lepton sector for the case of a Dirac neutrino. 

The utilization of the seesaw mechanism also fails to produce tribimaximal mixing. Since 
all the right handed neutrinos are all localized at the same place, all the overlap integrals 
which contribute to the Majorana mass matrix are of the same order of magnitude unless 
there is a substantial hierarchy amongst their domain wall parameters. Hence, for most 
solutions of interest, the right handed Majorana mass matrix assumes a non-hierarchal 
form. However, the neutrino Dirac mass matrices will maintain their hierarchical form since 
the lepton doublets will still be separated, and therefore the effective left handed neutrino 
Majorana mass matrix, ~ —m^M^mo is rendered hierarchical, and thus small lepton 
mixing angles for a Majorana neutrino will result. 

There are several approaches one could take to the problem of lepton mixing in this 
model. The most obvious is the inclusion of a discrete flavor symmetry like A4, or its double 
cover T'. This has in fact been employed successfully in RSI |10] and orbifold models |¥T| . 
Such an inclusion is beyond the scope of this paper. 

It was also not surprising that a solution was not found in the Dirac case, since the 
initial assumption that the couplings to 77 and x commuted cut the number of background 
Yukawa couplings to 15, and the assumption of universal electroweak Yukawa couplings 
cut the number of free parameters in the electroweak Yukawa sector to 1, giving 16 free 
parameters in total which determine the masses. If we want to generate everything except 
CP violation from non-hierarchical electroweak Yukawa couplings then, we need to generate 
the correct mass ratios for the charged fermions, quark and lepton mixing angles, the correct 
Am^2 ^^"^ -^^23 ^^"^ acceptable neutrino mass scale, which amounts to 18 constraints. 
Thus, we were never guaranteed such a solution. Letting go of the initial assumption that 
the background couplings commute allows us to introduce mixing angles and CP phases 
from that sector. In practice, solving the relevant equations for the fermion profiles with 
non-commuting and /i^^, is difficult and must be solved numerically; such an analysis 
will be deferred to a later paper. 
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IV. CONCLUSION 



As we have seen, the utihzation of extra dimensions is particularly useful in explaining the 
fermion mass hierarchy problem. This is particularly evident in the analysis without quark 
and lepton mixing, where we were able to show that the mass hierarchy which spans at least 
14 orders of magnitude could be generated from a set of domain wall Yukawa parameters 
which have a spread of roughly an order of magnitude. Furthermore, this spread could be 
reduced even further by making the Higgs profile more localized. As an added bonus in this 
analysis, by choosing parameters such that the colored Higgs was well displaced from the 
domain wall and the electroweak Higgs, the doublet-triplet splitting problem was solved and 
proton decay suppressed to such an extent that the colored Higgs mass no longer had to be 
fine tuned. 

Generation of quark mixing from the overlaps after initially assuming non-hierarchal 
SU (5) electroweak Yukawa coupling constants also looks promising, and we successfully 
generated the Cabibbo angle and the fermion mass hierarchy for the case with two gener- 
ations and Majorana neutrinos. Generically, small quark mixing angles and fermion mass 
hierarchies naturally arise from hierarchical mass matrices, although a more thorough nu- 
merical analysis will have to be done to find solutions for the full CKM matrix. 

We have given some arguments as to why the problem of tribimaximal lepton mixing 
problem cannot be solved simultaneously with the quark mixing and fermion mass hierarchy 
problems in this braneworld model. Typically with solutions to the latter two problems, the 
lepton doublets are spread out away from the Higgs profile, rendering the electron and Dirac 
neutrino mass matrices hierarchical, leading to small mixing angles. We believe this may 
be amended with the addition of a discrete fiavor symmetry like /I4, or by dropping the 
assumption that the 7] and x couplings commute. 

The addition of a flavor symmetry to the model as well as a more thorough analysis of 
the parameter space will be treated in later papers. We also must make a number of other 
additions to the analysis of this paper, most notably the inclusion of gravity and an analysis 
of the renormalization group evolution of the mass parameters. 
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